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We examine the particle production during tunneling in quantum cosmology. We 
consider a minisuperspace model with a massive, conformally coupled scalar field and 
a uniform radiation background. In this model, we construct a semiclassical wave 
function describing a small recollapsing universe and a nucleated inflating universe 
("tunneling from something"). We find that the quantum states of the scalar field 
in both the initial and the nucleated universe are close to the adiabatic vacuum, 
the number of created particles is small, and their backreaction on the metric is 
negligible. We show that the use of the semiclassical approximation is justified for 
this wave function. Our results imply that the creation of the universe from nothing 
can be understood as a limit of tunneling from a small recollapsing universe. 



I. INTRODUCTION 

A semiclassical picture of quantum cosmology based on the Wheeler-DeWitt equation 
describes tunneling from a state of vanishing size ("tunneling from nothing") to a closed 
inflating universe. The process of tunneling from nothing can be thought of as a limit of 
tunneling from a closed recollapsing universe of very small but nonzero size to an inflating 
universe ("tunneling from something"). This paper continues the investigation of particle 
creation during tunneling in quantum cosmology. Conflicting claims of excessive particle 
production that invalidates the semiclassical approximation Jl|,l2| , on the one hand, and of 
essentially no particle content in the nucleated universe [3, y,^, on the other hand, have 
been advanced in the literature. Our intent is to resolve this long-standing controversy. 

The process of tunneling from a recollapsing universe sensitively depends on the quantum 



state of that universe. In the companion paper 



we have shown that, at least in the case of 



a massless field, the results of Rubakov et al. Q,|2|| should be interpreted not as an indication 
of a large particle production but as a consequence of an inadequate choice of the initial 
quantum state of the universe. A generic quantum state of the recollapsing universe will 
contain a superposition of various semiclassical geometries. Some geometries in the super- 
position do not describe a nucleated universe but rather a universe that expanded from zero 
size "over the barrier" without tunneling, because it contained a large number of particles. 
Other geometries in the superposition will describe a nucleated universe with a small number 
of particles. One can hope to extract well-defined particle numbers from a Wheeler-DeWitt 
wave function only if it describes a single semiclassical background spacetime. However, 
the process of tunneling amplifies the differences between branches of the wave function. A 
semiclassical branch of the wave function describing a universe with a large particle content 
may be strongly suppressed in the initial recollapsing universe but may dominate the wave 
function of the tunneling universe (because it avoids being exponentially suppressed during 
tunneling). 

For a massless field, there is a well-defined vacuum and the particle production is absent 
0| . The quantum state chosen by Rubakov [jj corresponds to a squeezed state with respect 
to that vacuum and consists of an infinite superposition of excited states with all possible 
energies. Each excited state is described by a branch of the wave function that has its 
own semiclassical geometry. As we showed in Ref. 0|, one cannot neglect the difference 
between these semiclassical geometries if one considers the nucleated universe. We call this 
phenomenon "critical branching". We have shown that the "catastrophic particle creation" 
found by Rubakov et al. is in fact a sign of critical branching that happens with their choice 
of the initial state. If one chooses the initial state to be the vacuum, or any excited state 
that is an eigenstate of the particle occupation numbers, then there is no branching, and 
the wave function describes a single semiclassical geometry. 

n 

The purpose of the present paper is to extend the analysis of Ref. [6] to the case of a 
massive field. A nonzero mass m introduces important qualitative differences into the prob- 
lem. First, a massive scalar field is nontrivially coupled to the time-dependent metric, and 
there is necessarily some particle production. Even if one imagines starting the recollapsing 
universe in a vacuum state at early times, the field will be in a superposition of states with 
all possible occupation numbers at other times. Additional creation of excited states may 



occur as a result of tunneling. Thus, it is not clear a priori that critical branching can be 
avoided for any choice of the initial state. Moreover, the definition of the vacuum state for 
a massive field in an expanding universe is notoriously ambiguous; as a result, the number 
of particles is also subject to ambiguity. In this paper we give a detailed analysis of these 
issues. 

The approach taken by Rubakov et al. was to define the vacuum state by diag- 

onalizing the scalar field Hamiltonian at a moment of time. However, this procedure is 
problematic since it is known to give unphysical results for the particle density in some 
cases A well-motivated definition of vacuum in a slowly changing background geometry 
has been developed by Parker and Fulling; it is the so-called adiabatic vacuum p, 0|. An 
adiabatic vacuum of A;-th order is defined using a truncated asymptotic series obtained from 
the WKB expansion. It has some dependence on the moment of time when it is defined and 
on the order k at which the series is truncated. The resulting uncertainty is of order of the 
k-th power of the adiabatic parameter. On the other hand, it follows from our analysis in 
[f| that an admissible state of the universe must be defined with an exponential accuracy 
to avoid the "critical branching". We need to specify a vacuum state with a higher accuracy 
than the definition of the adiabatic vacuum allows. 

Our main result in this paper is a prescription for constructing a well-behaved semiclas- 
sical state of the universe which does not exhibit critical branching. We show that this state 
always exists and coincides with the adiabatic vacuum within the accuracy to which the 
latter can be defined. It also coincides with the canonically defined vacuum in the massless 
case. 

The paper is organized as follows. In Sec. HU we construct a family of well-behaved 
semiclassical wave functions of the universe. We give a prescription to select a unique 
quantum state that we call the "Gaussian vacuum". In this state the backreaction of produced 
particles is negligible and there is no branching. We obtain approximate expressions for the 
wave function of this state. In Sec. Mil we give a particle interpretation of the Gaussian 
vacuum state. We show that the Gaussian vacuum is indistinguishable from an adiabatic 
vacuum state to all allowed orders of the adiabatic expansion. In Sec. HXl we discuss our 
results as well as the issues raised by Rubakov et al. We also compare our results to those 
of Bouhmadi-Lopez et al. [1Q| who addressed some of the same issues. In the Appendices 
we give technical details of the calculations and, in particular, check the validity of our 



approximations. 



II. THE GAUSSIAN VACUUM STATE 



We consider a homogeneous closed FRW metric and an inhomogeneous massive confor- 
mally coupled scalar field. The classical action of the system is 
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Here R is the scalar curvature and the parameter H represents the cosmological constant. 
We use the Planck units, G — h — c — 1. We assume ff < 1 and m < 1 as a typical 
physically motivated case. 

We expand the inhomogeneous scalar field in spherical harmonics on the 3-sphere, 

7T\/2 



a(t) 



(2) 



n,l,p 



A rescaling by the factor a(t) is done for convenience. Below, only the index n will enter 

the equations, and we suppress the indices I, p of the modes Xnipif)- The summation over 

degenerate indices /, p spans / = 0, n— 1 and p = —I, / and introduces an extra factor 

n 2 which we shall insert in explicit calculations below. 

The wave function of the universe depends on all modes Xn of the scalar field, ^ = 

^ ( a ) {Xn})- The Wheeler-DeWitt equation, after appropriate rescalings Q|, takes the form 
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Here and below we explicitly write the Planck constant h = 1 only as a formal bookkeeping 
parameter, to clarify the use of the WKB approximation. 

In addition to the scalar field we now include a small amount of homogeneous radiation 
with energy density 

p r = a~ 4 e r , (4) 
where e r > is a constant parameter. Then the Wheeler-DeWitt equation is modified, 
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where we have defined 



V{a) = a 2 - H 2 a\ (6) 
cu n (a) = Vn 2 + m?a 2 . (7) 

If restricted to the coordinate a, Eq. Q is a stationary Schrodinger equation for a particle 
in a potential. 

A. Gaussian solutions of the WDW equation 

In the companion paper [f| we have used the method of Refs. [3, Q| to find 
an approximate solution of Eq. © . The solution may be found as a linear combination of 
Gaussian terms of the form 



* 0, {Xn}) = exp 
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where S(a) and S n (a) are functions to be determined. The functions S n (a) must satisfy the 
condition that we shall call the "regularity condition", 

< Re S n (a) < +oo for all a, n. (9) 

With this condition, the Gaussian wave function of Eq. ® is well-defined everywhere and 
quickly decays at large Xn- In Ref- [6] we have shown that a violation of the regularity condi- 
tion indicates a splitting of the wave function into decoherent branches with macroscopically 
different semiclassical geometries and different particle contents (critical branching). A wave 
function of the universe can be consistently interpreted in terms of a classical spacetime with 
a quantum field only if a single underlying semiclassical geometry is present. Each term of 
the form ((HI) will describe a single semiclassical geometry if the regularity condition holds for 
its function S n (a). Then the branches of the wave function corresponding to each such term 
will describe independent, decoherent semiclassical universes. Therefore we may impose the 
regularity condition on all terms of the form of Eq. (JBJ) that comprise the wave function of 
the universe. 



We now substitute the ansatz of Eq. (JHJ into Eq. (JHJ and obtain 

{S') 2 - V(a) - hS" + h S n = 0, (10) 
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Our approximation consists of disregarding the terms of order 0(h) and 0(%£); the ap- 
plicability of this approximation is analyzed in Appendix |Bj We then have the following 
equations for the functions S(a) and S n (a), 

(S') 2 -V(a)+e r =0, (12) 
S'S' n -Sl + ul =0. (13) 

It is clear that exp(— S/h) is a standard WKB wave function for a particle in a potential 
V(a) with energy e r . For H 2 e r <C 1 the turning points ai )2 are approximated by 



a? = 2H2 « en (14) 
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If 4H 2 e r < 1 then ai ^ a 2 are real and there exist two classically allowed regions and a 
classically forbidden region (see Fig. P). The physical picture of the resulting cosmology is 
the following. A small closed universe filled with radiation of energy density e r a~ 4 expands 
until the maximum scale factor a\ and then recollapses to a = 0. An expanding universe 
is created with a scale factor a 2 by tunneling through the potential barrier. If the WKB 
approximation is valid, then the scale factor a is a semiclassical variable and any monotonic 
function of a can be used as a time coordinate in the two classically allowed Lorentzian 
regions. The turning point a = a 2 corresponds to the beginning of time in the nucleated 
expanding universe; the Euclidean region a± < a < a 2 does not correspond to a classical 
universe. If H = 0, the second turning point is absent (formally a 2 — > +oo) and there is 
only one Lorentzian region < a < a\. 
We may rewrite Eqs. (|T2l - ([T31) as 



S(a) = ± J ^/V(a)-e r da, (16) 
±^V(a)-e r S' n = SI-loI (17) 
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Figure 1: "Tunneling from something": creation of an expanding universe by tunneling. 



In the two Lorentzian regions [0 < a < a\ and a > 02] the square root in Eq. fTKj) as- 
sumes imaginary values (with a positive imaginary part) and the upper sign denotes a wave 
exp[— S(a)) traveling to the right, with 

S(a) = iJ v^r - V(a) da. (18) 

In the Euclidean region \a\ < a < CZ2] the square root in Eq. (jTK|) is real and the upper sign 
denotes an exponentially decaying mode exp[— S(a)} with 

S(a) = J y/V(a) -e r da. (19) 

We shall denote by S n (a) and S~(a) the solutions of Eq. (fTTj) with the upper and the 
lower signs respectively, omitting the + subscript in S£(a) for brevity. If S(a) is given by 
Eqs. (jT%jl - (jTT?j) . then a linear combination 
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is an approximate solution of the Wheeler-DeWitt equation. The coefficients C± and the 
solutions S n , S~ must be chosen separately in each of the three physical domains (two 
Lorentzian regions and one Euclidean region). To match these solutions across the turn- 
ing points, we need to use the matching conditions of Ref. These conditions state, 
in particular, that the value of S n must be continuous across the turning points and 
moreover 

(21) 
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A similar matching condition must be satisfied at the second turning point a 2 , 

lim S n (a) = lim S~(a) = lim S n (a). (22) 

a— >a2— a^a2— a— ><J2+0 

Due to the tunneling boundary condition at a — > +00 (no wave traveling to the left) the 
coefficient C_ for the region a > a 2 must vanish and we do not need the branch S~(a) in 
that region. In addition to Eqs. (j2H)-(J221) we require that the regularity condition [Eq. (jHJ)] 
holds for S n (a), S~(a). We shall refer to such solutions of Eq. (fTTj) as "regular". 

Note that the function S~(a) does not have to be continuous at a — a\. This is because 
the branch with the lower signs [the second term in Eq. (f20|) ] is exponentially small at a = ai 
compared with the other branch and cannot be kept in the WKB approximation. 

B. Construction of regular Gaussian solutions 

Our task here is to obtain an explicit solution of 

VWTTr^ = Sl- ul (23) 
subject to the regularity condition [Eq. (jOJ)], 

< Re S n (a) < 00 for all a. (24) 

Here S n (a) is in general a complex function. Because of matching conditions, S n (a) must 
be continuous across the turning points a = a% t 2- Equation (|23|) is identical to Eq. (|TT|) : the 
equation for S~(a) can be obtained by changing the sign of the derivative d/da in Eq. (|23|) . 
Following Refs. [ill, Q, Ql, we introduce the semic assical time variables: t in the clas- 



sically allowed regions and r in the under-barrier region, 



— = \J e T — V(a), a < a\ or a > 02, (25) 
da 



— = ^/V(a) — e r , ai < a < a 2 . (26) 
dr 

The variable t is the conformal time in a Lorentzian universe, while r is the Euclidean 

conformal time. Below we shall be using the variables a, t and r interchangeably in the 
appropriate ranges of a. Then Eq. (|23|) becomes 

.dS n 2 2 „ ^ /ota 

1— — = b n — uJ n , a < ai or a > 02, (27) 
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= b n — u n , ai < a < a 2 . (28) 

dr 



The equations for the other branch S~ of Eq. (|T7jl are obtained by reversing the signs of 
the time derivatives. The branch S n corresponds to the wave traveling to the right in the 
classically allowed regions and to the decaying underbarrier branch in the Euclidean region. 

To analyze the existence of solutions and to obtain approximations, it is convenient to 
transform the function S n (a) and introduce an auxiliary function 

Ua) = (29) 

O n -f- LU n 

We imply that the pair of functions S^(a) is transformed into a pair Cn ( a ) bvrt will often 
omit the superscripts ± for brevity. The simple transformation of Eq. (J2*9"J) allows to use 
the functions and Cn interchangeably. The function ( n (t) can be interpreted as the 
"instantaneous squeezing parameter" [f| describing a state of the mode Xn that is related to 
the vacuum of the instantaneous diagonalization picture by a time-dependent Bogolyubov 
transformation. However, for our purposes it is enough to regard Eq. (j2~H| as a formal 
transformation that simplifies calculations. 

The regularity condition [Eq. (j^J)] is now equivalent to 

IC«(a)l < 1- (30) 

This condition is violated by |C n | = 1 when ReS n (a) = or by ( n = 1 when (S^a)! — > oo. 
From Eq. (|29|) and the matching conditions for S n (a) it follows that the (complex) function 
Cn(o) must be continuous across the turning points a 12 . From Eqs. (|2?j) . (|28|) we obtain the 
equations for ( n (a) in each region in terms of the conformal time variables, 

s«- = 2! - c " " i^r (1 " cS - (31 » 

Tr<^^-^ (1 -°- <32) 
Our goal now is to obtain a solution for ( n (a) satisfying the regularity condition. 

In AppendixElwe derive some technical results concerning regular solutions of Eqs. pTj) ~ 
Ip2]l. In Appendix IA H it is shown that a solution Cn( a ) is regular if |Cn(°2)| < 1, while the 
solution Cn ( a ) on ly needs to be regular at a — ax to be regular everywhere. From this it 
follows that regular solutions of Eq. (fTTjl always exist. For instance, one could numerically 
integrate Eq. (|17|) with the upper sign, starting at a = a 2 with any value S n (d2) such that 
< Re S n (a 2 ) < +oo, and obtain a regular solution at a > a 2 and at a < a 2 . 



To build a perturbative expansion, Eqs. (j31f - (|32j) can be rewritten as integral equations. 
For instance, the function („(t) in the Euclidean region, with an arbitrary boundary value 
(71(72), satisfies 

f>T 2 
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2ui r , 



uj^dr 
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(33) 



This equation can be solved iteratively starting with ( n (r) = 0. A similar calculation is 
performed for ( n (t) in the Lorentzian regions. In Appendix IA 21 we prove that the sequence 
of iterations always converges. 

An important case is an adiabatic (slow) expansion of the universe. The adiabaticity 
condition is 



1 dllJr. 



ui 2 dt 



m 2 a^/\e r — V(a)\ 



< 1. 



(34) 



(n 2 + m 2 a 2 ) 3 / 2 

Analyzing Eq. (|3*4l) with V(a) given by Eq. (JHJ), we find that it holds if n ^> m/H or if 
m ^> H. The only case when the ratio of Eq. (|3*4T) is of order 1 is when n ~ 1 and m ~ H . 

In Appendix IA 21 we also show under assumption of adiabaticity that there exist solu- 
tions for which |Cn(°)| i s always small, of order \u\j n /oj 2 n \ 1, and obtain the following 
approximations to general solutions [see also Eq. l)A14j) ]. 
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Cn( T ) = Cn( r l)exp 
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Here the given boundary values Cn( r 2) and (~ (r x ) must be sufficiently small, but are oth- 
erwise arbitrary. Similar expressions hold for the functions ( n (t) and (~ (t) in Lorentzian 
regions. 

Using the above approximations for ( n and (~ , we find S 1 ^ and find the Gaussian wave 
function of Eq. (f20|) . Calculations verifying the validity of the Gaussian approximation in the 
neighborhood of Xn = and the validity of the WKB approximation for a regular Gaussian 



solution are given in Appendix |Bj We conclude that the regular Gaussian wave function is 
a valid approximation to an exact solution of the Wheeler-DeWitt equation. 



C. Non-uniqueness of regular wave functions 

The solutions ( n , (~ of Eqs. IpjoTl - ljHoT) depend on the arbitrary boundary values Cn(7"2) 
and Cn ( r i) respectively. The choice of these boundary values is a priori constrained only by 
regularity and by matching conditions. The matching condition Cnfa) = C~ ( r 2) leaves one 
free parameter in the resulting wave function. 

However, this remaining freedom of choosing a regular solution of Eq. (fTTjl does not lead 
to an appreciable variation in the resulting functions S^(a) in the Lorentzian regions. The 
properties of Eq. (|T7jl in the Euclidean region a± < a < a 2 are such that, as long as the 
regularity condition holds, S~(a,2) is almost insensitive to the initial value S~(ai), and in 
turn, the value S n (ai) is almost insensitive to S n (a,2)- More precisely, the initial value S~(ai) 
may vary throughout a wide range but the final value S~(a 2 ) is constrained to be within an 
exponentially small interval. 

This can be seen directly from Eqs. (jHoTl - ljHoT) . For instance, the influence of the boundary 
condition Cn( r 2) on the value Cn( r i) is suppressed by an exponentially small factor 

= exp [-29 b ] . (37) 

Here 6 b is in the notation of Appendix IA 21 [see Eq. (|A5[) ]. The quantity 6 b is never small 
since a ~ H' 1 , V(a) ~ H~ 2 and 

T 2 uj n (a)da I m 2 

6b= Ljvw^rrr + ip- (38) 

We find that ^6^1 when n > 1 or m > . 

Since S n (a2) = S~(a,2), the value S n (ai) is also constrained to a narrow interval. Small 
changes of boundary values of the functions S n (a), S~(a) at the turning points a\p lead to 
small changes of their values in the Lorentzian regions. Therefore, different choices of the free 
parameter S~(ai) yield only exponentially small changes in the resulting solutions within 
the Lorentzian universes. The width of the possible range of solutions within a Lorentzian 
region is 
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A physically interesting simple case is that of a single closed universe, H = 0. In that 
case there is only one Lorentzian region < a < a\ and the wave function in the Euclidean 
region has only one branch, S n (r). In the limit H — ► 0, we find 9b — > +oo. Therefore the 
influence of the boundary condition S n (a,2) at 02 — > +00 is completely suppressed and the 
value S n (ai) is uniquely determined by regularity. Similarly, the regular wave function is 
unique in the limiting case e r — > 0, also because 9b — > +00 [the integral in Eq. (|38|) diverges 
near a = 0]. [This conclusion does not depend on the adiabaticity condition implicit in 
Eq. J2|; in Appendix El we prove the uniqueness of the regular solution under much 
weaker assumptions.] However, in the general case (H 7^ and e r 7^ 0) the wave function is 
not uniquely fixed by the regularity conditions alone. 

D. Prescription for a unique Gaussian vacuum 

In the previous section we have found a one-parametric family of well-behaved Gaussian 
wave functions that can serve as vacuum states. We may impose an additional constraint 
on the wave function to specify it uniquely. 

The extra condition is that the branch S~ (a) must be continuous across the first turning 
point a\. Together with other matching conditions this gives 

3- {ax) = S n ( ai ). (40) 

In Appendix IA 51 we show that this condition is always satisfied by a unique pair of regular 
solutions S n (a), S~(a). This proves the existence and uniqueness of the selected state. We 
shall call this state a "Gaussian vacuum". 

We view the additional constraint of Eq. (J4*0"jl as a prescription for defining a unique 
regular vacuum state. The wave function resultingfrom this prescription has the advantage 
that it agrees with the vacuum obtained in Ref. [6] for the massless scalar field. All other 
regular solutions obtained in Sec. Ill CI differ from the Gaussian vacuum by an exponentially 
small amount (in Lorentzian regions). 

Note that the Gaussian vacuum prescription is not local, since it requires knowledge of 
the behavior of the of the frequency oo n (a) and of the potential V(a) everywhere under the 
barrier. 



III. INTERPRETATION OF THE GAUSSIAN VACUUM AND PARTICLE 

PRODUCTION 



A semiclassical Gaussian wave function is interpreted by making a transition from the 



Schrodinger picture of minisuperspace to a QFT in curved spacetime [Hi 112, H3| • Nat- 
urally, this can be done only in a Lorentzian region where a semiclassical spacetime exists. 
One introduces the conformal time t via Eq. (|25|) . The scalar field x is canonically quantized 
using the standard creation and annihilation operators a n , a^ for the modes \ n [here the 
extra indices I, p of the modes are again suppressed for brevity, cf. Eq. (j2J)]. The ampli- 
tude Xn of the n-th mode is promoted to an operator Xn and decomposed into creation and 
annihilation operators, 

Xn(t) = BL n u n (t)+aii;*(t). (41) 

A vacuum state is defined by a particular choice of the mode functions z/ n (t). If the 
Gaussian wave function is given by Eq. ((HI) with a certain set of S n (a) for all n, then the 
corresponding mode functions v n are found from 

S n (t) = ——f--, < a < ai or a > a 2 . (42) 

The interpretation of the Gaussian wave function [Eq. ©] in the Lorentzian regions is that 
the quantized scalar field is in the vacuum state defined by the above mode functions v n (t). 

The mode functions v n are determined by Eq. fi2f up to an arbitrary constant factor. It 
follows from Eq. (|27j) that v n (t) satisfy 

+ = 0. (43) 

This is the usual equation for the mode functions of the n-th mode of a (rescaled) conformally 
coupled scalar field. The mode functions may be normalized by the Wronskian condition 

i. (44) 



A. The Gaussian state and the adiabatic vacuum 

In the companion paper we have shown that a Gaussian wave function satisfying the 
regularity condition can be interpreted as a single semiclassical spacetime with a quantized 
field in a certain vacuum state. In Sec. Ill Dl we have defined a unique vacuum state (the 



"Gaussian vacuum") that satisfies the regularity condition. Now we need to compare this 
Gaussian vacuum with a physically motivated vacuum in our Lorentzian universes. 

Neither the spacetime of the recollapsing universe nor that of the nucleated universe 
possess a well-defined asymptotically static regime. This does not permit one to define 
unique "in" or "out" vacuum states of the QFT in these spacetimes. However, one can define 
an approximate "adiabatic vacuum" of the QFT in a nonstationary FRW spacetime 0| if 
the adiabatic approximation is valid for the mode functions of the field. We now show that 
the adiabatic vacuum in fact coincides with the Gaussian vacuum state we defined (up to 
a small correction). The difference between the two vacua is in any case smaller than the 
uncertainty inherent in the definition of the adiabatic vacuum. 

We assume that the adiabaticity condition holds; its precise formulation is 

dUn « 1 (45) 



u 2 



dt 

[see also Eq. (|A6J) ]. In that case one can use at least a few terms of the WKB expansion for 
Eq. To obtain the WKB expansion, one can introduce a formal parameter A, define 

uj n {t) = uj n (Xt) and use the ansatz 



W{t)dt 



(46) 



The WKB function W(t) is found as an asymptotic series in A . The first few terms are 



v 4<4 8o£, 

Each derivative of u n adds a power of A; at the end we put A = 1. The true small parameter 
of the expansion is the slowness of change of ui n (t), formalized through \u n \ <C uj\ and 
analogous conditions on higher time derivatives. 

The definition of an adiabatic vacuum depends on the chosen adiabatic order k and on an 
arbitrary fiducial time to. Let the function W^ k '(t) give the WKB series [Eq. (137}) ] truncated 
up to terms O (A 2fc ). The mode function u n (t) describing the adiabatic vacuum are required 
to coincide with the WKB solution of order k at t = t . The condition for this is 



± + iW (k) + ld lnWik) 
Vr, 2 at 



= 0. (48) 

t=to 



In an adiabatically changing background, the adiabatic vacuum of order k gives an apparent 
particle creation rate of order k + 1 in the adiabatic parameter. For a given metric the WKB 



solution is usable only until a certain finite order k max after which the WKB series starts to 
diverge. 

The condition of Eq. (|48l) fixes the value S n (t ). [For comparison, the instantaneous 
diagonalization approach sets S n (t ) = u n (to).] In Appendix IA3l we show that the asymp- 
totic series obtained for S n (t) in the WKB approximation is the same (to all orders) as the 
asymptotic series for S n (t) obtained from the Gaussian vacuum. Therefore, the Gaussian 
vacuum coincides with the adiabatic vacuum within the accuracy of its definition. 

This coincidence is not specific to the Gaussian vacuum. Any other regular Gaussian wave 
function as found in Sec. Ill CI will give a value S n (to) different by exponentially small terms 
of order exp(— 29b). The WKB series contains terms of order [ujn/uj^ and is insensitive to 
exponentially small contributions. In other words, the definition of an adiabatic vacuum of 
order k at time t contains inherent uncertainties of much larger magnitude. We conclude 
that the regularity of the wave function specifies a quantum state of the universe that is 
indistinguishable from an adiabatic vacuum of any applicable order. 

B. Quantum state of the nucleated universe 

Previous work suggests that in the case of tunneling from nothing (e r = 0) the scalar field 
in the nucleated universe should be in the Bunch-Davies (BD) vacuum state. Now we can 
consider the case of tunneling from something (e r ^ 0) and compare the Gaussian vacuum 
state in the asymptotic region a — > +oo with the BD vacuum. 

In the Schrodinger picture of minisuperspace, the mode functions u n (t) of a vacuum 
are determined by the solutions S n (t) [Eq. (|4*2j)]. Suppose that another vacuum state is 
determined by another set of solutions S n (t) and the corresponding mode functions u n {t)- 
The two vacua are related by a Bogolyubov transformation, 

v n {t) = av n {t)+l3v* n {t). (49) 

We can express the Bogolyubov coefficients a n , (3 n directly through the functions S n (t) and 
S n (t), as follows. If both sets of mode functions are normalized [Eq. (jSJ)], then ot n and j3 n 
satisfy \a n \ 2 — \/3 n \ 2 = 1. We can select any value oft and express u n (t) and u n (t) according 



to Eq. i|3njl . The solution is 



«n = *>>n - U*U n , (50) 
/?n = -*>rA + V n V n . (51) 

Then the physically measurable quantity \/3 n \ 2 that gives the mean occupation number in 
the mode Xn is found (again, independently of a fixed value of t) as 

s n {t) - s n (t) 2 

l&l = a 7^ „ ^ - (52) 



4i?e S^ite S n (t) 

In the case e r = the regular solutions S^(a) of Eq. (JT7l) are unique and are known to 
correspond to the de Sitter-invariant BD vacuum jj, 0|. The difference between the cases 
e r = and e r ^ can be seen from Eqs. PHI) . (|3H|) : the neighborhood of the first turning 
point r = 7~i contributes an exponentially small amount to S~(t2). Therefore the regular 
solution in the e r 7^ case [denoted for now S n (a)] may differ from the BD solution S n (a) 
only by an exponentially small contribution of order exp(— 29 b ). From Eq. (|52|) we find that 
the Bogolyubov transformation relating the two vacua is exponentially close to the identity 
transformation (a n — 1, (3 n = 0). 

We conclude that the Gaussian vacuum state derived from the regularity of the wave 
function contains an exponentially small number of particles from the point of view of the 
BD vacuum. 



IV. DISCUSSION 



We have shown that a minisuperspace model of quantum cosmology with a conformally 
coupled massive scalar field admits a vacuum state described by a well-behaved Gaussian 
wave function (the "Gaussian vacuum") . Such a vacuum state is provided by a regular Gaus- 
sian solution of the form of Eq. (JHJ) and represents a recollapsing universe and an expanding 
universe that is nucleated by quantum tunneling. The obtained wave function describes 
a single semiclassical geometry throughout the Lorentzian and Euclidean regions and does 
not exhibit any "branching" into different semiclassical geometries. We have checked the 
consistency of the WKB approximation, of the Gaussian approximation, and of neglecting 
the backreaction of the scalar field perturbations. We have also shown that the Gaussian 



vacuum describes both the recollapsing and the expanding universe in quantum states that 
are indistinguishable from an adiabatic vacuum of any meaningful order. 

Unlike the definition of an adiabatic vacuum, our prescription for the Gaussian vacuum 
is not local. In other words, we cannot specify our preferred quantum state at some value of 
a < ax without a full knowledge of the shape of the barrier V(a) and of the function u n (a) 
under the barrier. Nevertheless we believe that our choice of the vacuum state is adequate, 
for the following reasons, (i) Our prescription is in the spirit of the definition of an adiabatic 
vacuum which was designed to minimize the apparent particle production, (ii) The Gaussian 
vacuum coincides with any adiabatic vacuum within the inherent uncertainty of the latter, 
(iii) In the massless case, where a well-motivated independent definition of the vacuum 
is available, our prescription selects the correct vacuum state, (iv) Finally, the Gaussian 
vacuum represents a fixed semiclassical background geometry. Any excited states built from 
the Gaussian vacuum using a finite number of creation operators are also well behaved and 
represent single geometries. If we wished to consider tunneling from a different state of the 
recollapsing universe, e.g., from a squeezed state, it would be more illuminating to represent 
that state as an infinite superposition of excited states built over the Gaussian vacuum, 
with each branch representing a single semiclassical geometry. The physical interpretation 
of the QFT in resulting spacetimes would be unambiguous. It would be interesting if the 
Gaussian vacuum (or some state exponentially close to it) could be specified by a set of local 
conditions, but at present we are unable to provide such a specification. 

Our construction of the Gaussian vacuum is general and not specific to the particular 
physical system we considered. The recent work of Bouhmadi-Lopez et al. 1^| discusses the 
quantum cosmology of a FRW universe filled with massless radiation field and a conformally 
coupled massive scalar field, in the presence of a positive or a negative cosmological constant. 
They have used the regularity conditions to constrain the Gaussian wave functions and 
obtained a family of regular Gaussian solutions. This family corresponds to our family of 
regular solutions S n (a), S~(a) parametrized by the boundary value S~(ai) [Sec. Ill Qj . We 
have shown in general that all such regular solutions describe practically the same state of 
the Lorentzian universes (up to exponentially small corrections). Our considerations are 
more general than those of Ref. [l^, where the existence of a family of regular solutions was 
demonstrated explicitly in a particular model with a negative cosmological constant. 

Our analysis can also be applied to artificial situations such as those treated numerically 



by Rubakov et al. in Ref. |2|. We may consider Eq. (J2*3*jl with arbitrary functions V(a), 
uj n (a) as long as the potential V(a) has the same qualitative behavior as the function of 
Eq. (JBJ), providing a potential barrier. Our statements about the existence and the behavior 
of regular solutions hold for any suitably well-behaved functions V(a) and u n (a). The unique 
Gaussian vacuum can be found either by a perturbative expansion or numerically. 

In Ref. jj| the functions w n (a) and V(a) were selected so that the resulting Lorentzian 
universe had an asymptotic region a — > — oo, where 0J n (a) — > uj^ = const. The asymptotic 
in- vacuum is then well-defined and is specified by S n (a — > — oo) = Un\ However, the 
Gaussian vacuum for such situations will generally specify S n (a) such that 

lim [S n (a) - u n (a)] = AS n ^ 0. (53) 

The quantity AS n can be found using the methods of Appendix IA 21 If the adiabaticity 
condition is satisfied, all AS n are exponentially small. (Note that the potentials used in 
Ref. P| do not satisfy the adiabaticity condition.) 

Nonzero values of AS n imply that the in-vacuum is an infinite superposition of semi- 
classical wave functions with different occupation numbers over the Gaussian vacuum. In 
fact, it is a squeezed state with squeezing parameters ( n ~ AS n /2uJn\ Even if states with 
large particle numbers are exponentially suppressed to the left of the barrier a < a\, critical 
branching may still occur. The under-barrier suppression of lower-energy branches is also 
exponential, and thus one is faced with a quantitative question of which of the two effects 
prevails. The answer to this question is likely to be model-dependent, and it is conceivable 
that a wave function starting in an in-vacuum state in the asymptotic region can exhibit 
critical branching under the barrier and be dominated by highly excited branches after the 
tunneling. We shall not attempt to address this issue in the present paper. 

In the physically realistic case, when there is no asymptotically static region, the Gaus- 
sian vacuum can serve as a suitable definition of the vacuum state. The wave function for 
"tunneling from nothing" can then be obtained as a limit of tunneling from a small recollaps- 
ing universe when the energy of that universe vanishes, e r — > 0. In this limit, the Gaussian 
vacuum becomes the de Sitter invariant Bunch-Davies vacuum state. 
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Appendix A: PROPERTIES OF A REGULAR VACUUM SOLUTION 



1. The regularity condition 

Here we analyze the regularity condition [Eqs. (JOJ and (|30|) ] and demonstrate that the 
regularity condition for all a is equivalent to imposing the regularity condition only at the 
turning point. First we show that if |Cn(^o)| < 1 at some t = t within a Lorentzian region, 
then |Cn(£)| < 1 for all other t. Then we demonstrate that the condition |Cn(V)| < 1 will 
hold in the whole Euclidean region a± < a < a 2 if it holds at the turning point a 2 . [For the 
other branch (~ (r), a similar argument will show that the regularity condition at a — a\ is 
sufficient.] No assumptions of adiabaticity of u n (a) are made. 

The function ( n (a) is a differentiable (complex) function satisfying a first-order differential 
equation with continuous coefficients, and a Cauchy problem will have a unique solution. 
Therefore each point £n in the two-dimensional configuration space (the complex ( n plane) 
has a unique solution ( n (a) starting from ffi at some a = a . 

First, consider a Lorentzian region (either < a < a x or a > a 2 ). The function ( n (t) 
satisfies Eq. (|3Tf and therefore 



CCn 



The quantity 2Re 



Re [CCn] = ~ ICnl 2 = [1 - |Cn| 2 ] Re C (Al) 

is the "radial velocity" at a point ( n . It follows that any trajectory 
Cn{t) starting on the circle \( n \ = 1 will remain on that circle. From uniqueness, it follows 
that no trajectory can cross the unit circle: any solution ( n (t) that is not entirely on the 
circle is either completely inside or completely outside of the circle. 

Second, consider a Euclidean region where ( n ( T ) satisfies Eq. (p2j) . A similar calculation 
for |Cn| = 1 gives 

I— \Cn\ 2 = 2uj n \Cn\ 2 > at iCnl = 1. (A2) 

2 GST 



Therefore all trajectories Cn( r ) that cross the unit circle must go outwards at the crossing 
point. The condition |Cn( r 2)| < 1 guarantees |Cn( r )| < 1 for all r < r 2 . If we take r 2 to be 
the value of r at the second turning point a = a 2 , then the required statement |Cn(o)| < 1 
for di < a < a 2 follows. 



2. An approximate solution 



The main statement of this section is the following. If the adiabaticity condition 



1 d,Ur. 



2u 2 dt 



< 1, 



1 dUr, 



2uj 2 dr 



<1, 



(A3) 



is satisfied and U(t), L0(t) are slowly-changing and non-oscillating functions, then there exists 
a regular solution ( n of Eqs. (|3*Tj) . ([3^1) which is always small, \( n \ <C 1. The proof is 
constructive and uses an integral equation for ( n to build a perturbative expansion. The 
method is based on a time-dependent perturbation theory; a somewhat similar but more 
cumbersome treatment is in Ref. As a by-product we shall find all regular solutions 
Cn(a)) parametrized by the boundary condition at a = a 2 . 

It will be convenient to consider one equation for ( n as a function of one complex variable 
instead of two Eqs. (pTTjl - ljH^l) using two different time variables. Start with Eq. and 
define the new "time" variable 9 by 



9(a) 



uj n {a) 



-da. 



(A4) 



'oa V £ r ~ V ( a ) 

Here the square root has the standard branch cut, Re y/z > 0, and the contour integration 
in Eq. (|A4J) uses a — > a + i8 with 5 > at the poles a = ai >2 . The value of 6 is real and 
positive for a > a 2 and coincides with the phase 9{t) = f* uo n dt. In the Euclidean region 
{a\ < a < a 2 ) we obtain, with the above branch cut, 9(r) = i J^ 2 u n dr. In the Lorentzian 
region < a < a,\ we have 9{t) = i6b + f t 1 oj n dt, where 

uj n (a)da 



LUndr 



y/Vic 



(A5) 



[a) — e r 

is the under-barrier "Euclidean phase" and ti >2 , ri j2 are the values corresponding to a = ai j2 . 
For our purposes it is enough to consider only the half-plane ImO > 0. 
Denote for convenience 



1 du n m 2 a\/e r — V(a) 
= 2(n 2 + m 2 a 2 fl 2 



2uj 2 n dt 



(A6) 



If the adiabaticity condition is satisfied, the value of / is everywhere small and we can find 
a bound / such that \f(a)\ < f <C 1. 
The equation for ( n (9) is 

^ = 2i( n -(l-£)f(9). (A7) 



From our analysis in Appendix IA II it follows that a solution ( n (a) which is regular, 
|Cn(o)| < 1 for all a, is uniquely determined by its value at the second turning point a = a 2 
(i.e. at 9 — 0). Therefore we assume a boundary condition ( n {6 = 0) = b with an arbitrary 
|fe| < 1. Then the solution of Eq. (|A7j) satisfies an integral equation, 

U9) = be™ - f e 2 ^- 9 ') (1 - Cm) f{0')d9'- (A8) 



o 

The integration in Eq. (|A8|) is understood as contour integration along e.g. a straight line 
connecting and 9. [The function /(a) has a branch point at a = in/m and we assume that 
an appropriate branch cut is imposed in the complex 9 plane.] 

We assume that \f(9)\ < f for the relevant values of 9, where fo is a fixed number. 
Then the integral equation for ( n (9) can be solved iteratively, starting with ( n = 0. This 
corresponds to a perturbation theory expansion in fo for Eq. (|A7j) . 

We now show that the iteration of Eq. (|A8|) always converges as long as Ira 9 > 0. Denote 
Cn\9) the k-th element of the iteration sequence. The initial function is 



C (°)(0 ) = be ™_ e H6-o') f (e>)d9\ (A9) 
Jo 

and the next approximations are found as 

e +1) (0) = (i \9) + f e^-e') [(i k \9')] 2 f(9')d9'- (MO) 
The difference between successive approximations is 

C ( fc+1 ) - C W= f dff f(9')e 2i ^ 



x([e ) (^)] 2 -K fc - i) (^)] 2 )- (ah) 

Now we shall estimate the integral in Eq. (|Allj) and show that the LHS tends to as k — > oo. 



From < Im6' < ImO we get |exp [2i{6' - 9)}\ < 1. Since \( n {9)\ < 1 and \f(6)\ < f , 
we can estimate 



< /, 



(i k+1 \9)-(i h \e)\ 



[Q 



(fc-1), 



d9' 



< 2/ / \^\e')-Ct l \o')\dd'. 



(A12) 



Starting from 



< /o |^|, we can prove by induction that 



\<£Xo)-<£-Ho)\<^^. 



2k\ 



(A13) 



This sequence clearly tends to as k — > oo at fixed 9. 

In terms of real time variables, Eq. (jA8jl can be rewritten as e.g. for the S~ branch in 
the Euclidean region, 



Cn 0) = Cn Ol) exp 



-2 



LU n dT 







/ exp 


-2 y o; n dr 







(A14) 



We have obtained the solution ( n (9) of Eq. (|A7ft for the boundary condition ( n (9 = 0) = b 
(with arbitrary |6| < 1) as the limit 

( n (9) = lim (W(9) (A15) 

of a sequence Cn fc) defined by Eqs. (fATHl- lfAToj) . 

Now we can show that if a boundary value is small, |6| fo, then ( n is always small 
and at most of order / in the adiabatic case /o< 1. The first iteration [Eq. (|A9|) ] clearly 



satisfies 



d°\9) 



fo; each subsequent iteration will only add terms of higher order in 
fo. Therefore, the dominant term of the solution is given by Eq. (fA^l) . Since Im9 > 0, 
it is clear from Eq. (|A9j) that |Cn(#)| is of order |/ | everywhere. We now only need to 
check that the oscillating integral in Eq. (|A9|) does not accumulate a large value when 9 is 
large. [This could happen, for example, in the case of a parametric resonance when oj{t) 
is oscillating.] By assumption f{9) is itself a slowly-changing function of 9 and we can 
approximate f{9) ~ /i + fi9 where f\ and fi = df jd9 are small constants of order f . Then 



the integral over one oscillation between 9 and 9 + 2n is 

2 \fi + h9)e- 2ie d9 = mf 2 . 



(A16) 



The value accumulated over many oscillations between 9\ and 9 2 = 6± + 2ixk can be approx- 
imated as 

X>|w + 2 ^Kf §«»= 5 l/«W (A17) 

This value is bounded by |/ |. Therefore the integral of Eq. (|A9|) remains small also for large 
values of 9. 



3. An asymptotic series for the adiabatic case 

We can use Eq. (|A8ft to obtain an asymptotic expansion for ( n (9) if the adiabaticity 
condition f <C 1 holds. We shall also assume that b ~ fo. The expansion is in the number 
of derivatives in 9 as well as in powers of / and is found through integration by parts, e.g. 



(2i) n + 1 dO n 



(A18) 



J* n=0 

Similarly to the argument leading to Eq. (|A12J) . one can show that the iteration sequence 
of Eq. (|A8|I gives a convergent sequence of asymptotic series in which the k-th iteration 
changes only terms of order k and higher. Therefore the asymptotic series for the solution 
Cn{9) is the limit of this sequence. [Each asymptotic series in the sequence, like the series 
in Eq. lfAT8|) . may not actually converge.] In the Lorentzian region < a < a\ we can omit 
the exponentially suppressed terms proportional to exp(2i#). Then the first terms of the 
resulting asymptotic series are 

(here the prime denotes d/d9 and all derivatives of / are evaluated at the same point 9). 
Converted back to the time variable t, this becomes 

t , A If u n u 2 n 

Qn[t) ~ — A — ~o \ 3 2 ~ 

16 \cj* u* ) [ 1 

This asymptotic series is a power series in the adiabatic parameter and necessarily misses 
any exponentially small contributions. 



4. Asymptotic series from the WKB approximation 



In the adiabatic case the WKB approximation may be applied to Eq. (|4*3|) to obtain a 
solution 



u n (t) oc 



1 



exp 



—% 



Wdt 



(A21) 



The auxiliary function W(t) is found as an asymptotic WKB series. One may iterate the 
equation 



W {k+i) 



1 WW 



\| ^ 2 WW ^ 

starting from W^ = co n (\t) and expanding in the formal adiabatic parameter A. At the end 
one sets A = 1. Then the solution S n (t) is obtained from Eqs. l|4*2*J) . (|4*6|l and transformed 
into Cn(^) using Eq. (f29|) . The first terms of the resulting series are 



WW 

ww 



(A22) 



S n (t) 
Cn(t) 



U n + 



2iuj 
i Co n 1 
4^ "8 





3^ 










^1 





+ ... 



(A23) 
(A24) 



It is clear that Eq. (IA24|) should coincide with the asymptotic series obtained above in 
Eq. (|A20J) using a different approach. An asymptotic expansion in powers of the adiabatic 
parameter will necessarily miss any exponentially small contributions, due to the nature of 
the power-law expansion. The remaining asymptotic series is unique for a given function 
(n{t), whether it was obtained from a WKB expansion or from any other procedure. An 
advantage of the method of Appendix IA 31 is that it can compute exponentially small terms 
in the solution. 



5. Existence and uniqueness of the Gaussian vacuum 

Here we prove that for any continuous function UJ n {r) on an interval [ti,t 2 ] there exists 
a unique pair of (complex) functions S n {r), S~(r) that satisfy the equations 

^ = % ~ (A25) 

= KT-<4 (A26) 
the regularity conditions for any r G [t 1 ,t 2 ], 

< Re S n (r) , Re S~ (r) < +oo (A27) 



and the matching conditions 



S-(r 1 ) = S n (r 1 ), S-(r 2 ) = S n (r 2 ). (A28) 

Consider an auxiliary function g(s,r) defined as the value g(s,r) = S n (r) obtained by 
solving Eq. (|A25jl with the boundary condition s = S n (r 2 ). Due to the uniqueness theorem 
for first-order differential equations, the function g(s, r) is differentiate and provides a one- 
to-one map of the complex s plane at any fixed r. It is also clear that g(s, r) has real values 
for real s. From our results in Appendix IA II it follows that < Reg(s,r) < +00 for any 
s such that < -Res < +00 and for t\ < r < r 2 . The final useful property of g(s,r) is 
\dg/ds\ < 1 for < Re s < +00. We can prove it as follows. The function dg(s,r)/ds as a 
function of r satisfies 

d dg dg dg 

TrYs =2g ^ T) ^ dS {T = T2) = 1 - (A29) 

This equation determines the function dg/ds at fixed s as 

dg, 



ds [T)=eXP 



T-2 



2 / g(s,r)dr 



(A30) 



Since Reg(s,r) > 0, we obtain \dg/ds\ < 1 for any r G [ti,t 2 ]. This means, in particular, 
that the map s — > g(s, Ti) decreases the distance between points in the complex s plane. 

Similarly, we define the function g~ by solving Eq. (|A26JI with the boundary condition 
S~(ti) = s to find S~(t) = g~(s,r). The function g~(s,r) has the same properties as 

g(s,r). 

The problem of finding S n , S~ is now equivalent to solving the algebraic equation 
g{g~(s,Ti),T 2 ) = s for s. The map s — > g{g~{s,Ti),r 2 ) clearly has the same proper- 
ties as the functions g and g~ . The existence of a real positive solution s > follows 
from the fact that the functions are continuous, have real values for real s, and satisfy 
< g(s,r), g~(s,r) < +00 for < s < +00. The solution is unique because if we assume 
that Si 7^ s 2 are two solutions, then the distance between Si and s 2 is decreased after ap- 
plying the map s — > g (g~(s, Ti), t 2 ), which is a contradiction since by assumption si >2 are 
stationary points of this map. [The unique solution can be obtained numerically by iterating 
the map.] 



6. Uniqueness of the regular solution for infinite barriers 



Here we consider the equation 

^ = S 2 -u 2 (r) (A31) 

GST 

on the interval < r < +00. We call a solution S(r) regular if < ReS(r) < +00 for 
all t > 0. The main statement is that Eq. (|A31|) has a unique regular solution S(t) if the 
function uj 2 {t) is continuous and bounded from below at large enough r, namely uj(t) > f 
at r > r , with a suitable constant / > 0. 

The technical conditions of this statement [continuity and a lower bound on u(r)] are 
sufficient but not necessary. The continuity of o>(t) is used only to ensure that the Cauchy 
problem for Eq. (|A31J) has a unique solution that is a continuous function of the initial 
condition. 

It is enough to consider the case co(r) > f for all r > because a unique regular solution 
in a domain r > r is uniquely extended to a regular solution for < r < t . 
Firstly, we prove that any regular solution S(r) must satisfy 

ReS(r)>f (A32) 

for all r > 0. For this it is enough to show that a regular solution satisfies ReS(r) > f — e 
for arbitrary e > 0. The differential equation for R(r) = ReS(r) is 

^ = R 2 - [Im S(r)} 2 - u 2 (t) < R 2 - f . (A33) 

So the function R{r) cannot grow faster than a solution Rq{t) of the equation 

^ = Rl ~ f- (A34) 

Given a regular solution S{r) and a point Ti, we can choose Ro(t) such that -Ro(ri) < 
ReS{ri) and then it follows that Rq{t) is a lower bound for ReS{r) for < r < n. Since 
we know that ReS{r) > 0, we can use the boundary condition i?o( r i) — at any T\ > to 
obtain lower bounds on ReS(r). The general solution of Eq. (|A34J) is 

1 - Ae 2fT 

R ^ = fTTA^F' (A35) 
where A is an integration constant. It is easy to see that any solution Ro(t) equal to zero at 
some large r must be exponentially close to / at smaller r. More precisely, for any e such 



that < e < f and for any T\ the solution Rq(t) will satisfy Rq(t) > f — e for < r < T\ if 
we impose the boundary condition Rq{t 2 ) = at 

r 2 = r 1 + i-ln^— ^. (A36) 

But the range of r is infinite and we can chose T\ to be arbitrarily large. Since Ro(r) is a 
lower bound for ReS(r), it follows that ReS(r) > f — e for < r < t\ with any T\ > 0. 
Therefore, we have shown that any regular solution of Eq. (|A31j) satisfies Eq. (|A32j) for all 

T. 

Secondly, we prove that the relative difference between any two regular solutions 5i(t), 
S2 (t) decreases exponentially with diminishing r. More precisely, for any two values r a , U 
such that r a < r bl 

\Si{r a ) - S 2 {r a )\ 2 < \Si(ri) - S 2 (T b )\ 2 ^_ 2f{Th _ Ta) 
IS^f + \S 2 (r a )\ 2 ~ IS.in)] 2 + \S 2 ( n )\ 2 ' 1 ' 

as long as these solutions satisfy Re Si t2 (r) > f for r a < r < r b . This can be proved directly 



by using Eqs. l|A31jl - (IA32l) to evaluate 



d ln IS^t) - S 2 (t)\ 2 



\Si\ 


2 ReS 2 +\ 


s 2 \ 


2 Re S 1 +u 2 Re(Si + S 2 ) 
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i(r 


)| 2 + 


S2(r)\ 2 



> 2f. (A38) 

From Eqs. (|A32j) and ()A37jl we find that the difference between any two regular solutions 
at r = must be equal to zero. This follows from Eq. (|A37J) with r a = 0: the relative 
difference by definition cannot exceed 1, and the RHS can be made arbitrarily small by 
choosing large enough r b . Therefore the regular solution is unique. 

The existence of a regular solution can be proved by construction. We find a regular 
solution S(t) as the limit of a sequence of solutions that are regular for a finite part of 
the interval < r < +00. Consider a (real) function S(t;t ) defined as the solution of 
Eq. (|A31|) with the boundary condition S (r ; r ) = 0. The function S (r; r ) is not a regular 
solution because it satisfies Re S(r) > only on the interval < r < r but not at larger r. 
We now show that the limit of S (r; r ) as r — ► +00 is a regular solution. The existence of 
the pointwise limit (taken separately at each r) 

S(t) = lim S (r; r ) (A39) 



follows from Eq. 1|A37J) : the relative difference between functions becomes exponentially 
small when r grows. It is clear that the resulting function S(r) will satisfy Eq. (|A32j) . 
It remains to show that the function S(r) defined by Eq. (IA39j) is actually a solution of 
Eq. (IA31jl . This follows from a continuity argument. A solution of Eq. (IA31jl is a continuous 
function of an initial condition. Therefore the limit of Eq. (JA39J) is the same as the solution 
of Eq. (|A31j) with the boundary condition 

S (r = 0) = lim S (r = 0; r ) . (A40) 

This completes the proof of the existence and uniqueness of the regular solution of Eq. (|A31J) . 



Appendix B: APPLICABILITY OF THE APPROXIMATIONS 

In this Appendix we analyze the applicability of the Gaussian and WKB approximations 
to the Wheeler-DeWitt equation. We use the WKB approximation in two places in our 
calculation. First, the WKB approximation is used in substituting the Gaussian ansatz into 
the Wheeler-DeWitt equation, when all terms of order 0(h) are ignored. Second, the WKB 
approximation is applied to Eq. (133*1) . We analyze the latter approximation first. 



1. Using the WKB approximation for Eq. (|43j) 

The WKB ansatz of Eq. (|A21j) is a good approximation for Eq. (f"""*"l) as long as 

w« < oj 2 n . (Bl) 

Comparing this with Eq. (|A20j) . we find that the squeezing parameter ( n (a) is small if and 
only if the WKB approximation is valid. 

Consider the behavior of ( n (a) in the Euclidean region where we use the time variable r 
(overdots will denote derivatives by r). Taking the leading term of Eq. (|A20J) . 



iti."ayJV(a) — e 



,2, 



Cn(fl) ^4u;2 - 4 ( n 2 +m 2 a 2)3/2 ( B2 ) 



and estimating V (a) — e < a 2 , we obtain 



m 2 a 2 



£™( a ) < T^^^^ — 2 2W2 - ( B3 ) 

4(n 2 + m 2 a 2 ) 6 / 2 ' 



At fixed n, this function reaches a maximum at ma = ny/2, which gives a bound 

C»(fl) < tA^- (B4) 
2ny 27 



Therefore |Cn( a )| <S 1 and the WKB approximation is valid if n ^> 1. Similarly, one can 
show that the WKB is valid if ma ^> n or if ma <C n. The only case when the WKB may 
not apply to Eq. (IB2ft is when n ~ 1 and ma ~ 1 at the same time; the maximum value of 
( n is at most ~ 0.1 in this case. 

In the Lorentzian region, Eq. ()B2|) gives at large a 

ICI ~ (B5) 

The adiabaticity condition is satisfied at large a only if H <C m. However, we know from 
Appendix IA II that the solution ( n (a) remains regular in Lorentzian regions even if the 
adiabadicity condition does not hold somewhere in those regions. 



2. Using the WKB approximation for the WDW equation 



When we substitute the Gaussian ansatz into the Wheeler-DeWitt equation [Eq. (jHJ)], 
we disregard terms of order 0(h). These terms are HSq, hS", and the "backreaction" term 
hJ2 n S n - The WKB approximation is valid if the disregarded terms are smaller than the 
typical magnitude of other terms in the respective equation. For the first two terms, the 
conditions are 

HS% < V{a), HS'; < col (B6) 
The WKB approximation may only be valid away from the turning points. We find 

I VI „o I1-2A 21 
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2^3/2 - " A 2(! _ A 2)3/2' 

this is small when H 1 away from turning points. We can also verify that 

2m 2 a~ 



(B7) 



Qll 



d 

da 



n 2 m 



3 



(B8) 



_uj n (a)_ 
for ma >n or ma n. 

Now we consider the backreaction of the excited field modes on the metric. The diver- 
gent sum hJ2 n S n corresponds to the infinite zero-point energy of the oscillators \n- We 
assume that this infinity is absorbed by an appropriate renormalization Q|. Any finite terms 



remaining after renormalization will be of order of the squared curvature ~ if 4 <C 1 and we 
neglect them here (see e.g. Ref. The remaining correction due to backreaction consists 

of replacing the effective potential 7(a) by 



7(a) + AV = V(a)-2hJ2 



n 3 p n , 



(B9) 



n=Q 



where p„ are occupation numbers in the modes Xn relative to an appropriate vacuum and 
we have inserted the multiplicity factor n 3 [see Ref. Jq], Eq. (37) for details]. The average 
occupation number in the state of the mode Xn characterized by Cn(a), relative to a vacuum 
described by (n\a), is found from Eq. ijH^J) . 



(Pn) = | A»| 





f AO) 




(1-Kn| 2 )(l- 


AO) 


2 

) 



(BIO) 



We shall not underestimate (p n ) here if we take (n^ = 0, which corresponds to the vacuum 
of the instantaneous diagonalization. The sum over large n can be estimated using the WKB 
approximation for ( n [cf. Eq. (IA20j) ] applicable at large n (as shown in Appendix IB lj) . The 
leading term is 

ICn| 



OJr. 



4^ 



(Bll) 

In this limit ( n is small and we can take 1 — \£ n \* w 1. At fixed a the squeezing parameter 
Cn(°) decays as n~ 3 at large n [Eq. (|B2|l]. This decay is fast enough so that the sum over 
modes in Eq. (|B9|1 converges. We obtain 

AV 



V 



7^1 

n=0 

oo 

E 



ic,i 



S 4 2 

n m a 



n=0 



n 2 + m 2 a 2 ) 3 



m 
32 



(B12) 



Since m < 1 in Planck units, we find that fractional change in 7(a) due to backreaction 
is small. [The occupation numbers are computed here in the instantaneous diagonalization 
vacuum where ( n is the instantaneous squeezing parameter. The particle numbers in a 
higher-order adiabatic vacuum are expected to be even smaller.] 



3. Validity of the Gaussian approximation 



In using the Gaussian ansatz, we disregard the terms which are quartic in \n but retain 
the terms quadratic in \n- The Gaussian ansatz gives rise to a quartic term 

v 4 / d <? \ 2 

Xn ( \ (B13) 



4 y da 

This term can be disregarded at small enough Xn if the following inequality holds, 

f(f) 2 «xM. (BM) 

The latter condition gives a corridor around the line Xn = in which the Gaussian ansatz is 
valid, \xn\ < Xmax- We can show that the width of this corridor is never small. The width 
of the corridor (divided by a, since Xn is rescaled by a) can be estimated using 

S n {a) ^uj n + ^- (B15) 

ZU! n 

and Eq. JA]|- We find 



X max \ 2 2u n (a) _ 2u n ^/V(a) 



n 



2 



dS n /da cl{S% — uj%) ~^ m 2 a 2 ^ 

This expression is bounded from below and therefore the allowed corridor never shrinks to 
zero width. 
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